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FORWARD 


The goal of Task 4/4a of NASA’s Large Engine Technology Contract (Contract NAS3-26618) 
has been to develop a design method for turbofan harmonic noise. The method is based on 
theory to the extent that it permits calculations at the high frequencies and wavenumbers 
required for modem engines. Grid requirements and high computer times preclude direct use of 
Computational Fluid Dynamic and, instead, techniques from classical duct acoustics and cascade 
unsteady aerodynamics have been adapted to the task. Based on earlier work, we felt that the 
new method had to include the effect of swirl between the rotor and stator and the effect of 
unsteady rotor/stator coupling. This lead us to a scheme in which the acoustic elements of a fan 
(inlet, rotor, stator, nozzle) of a fan are analyzed first as isolated elements and then coupled by 
matching eigenmodes at planes between the elements. The end product of this work is TFaNS - 
the Theoretical Fan Noise prediction and design System, which is reported separately by D. A. 
Topol: “Development and Evaluation of a Coupled Fan Noise Design System,” ALAA Paper No. 
97-1611. 

At the heart of the prediction system is the scattering behavior of the blade rows. This 
encompasses conversion of vortical waves (wakes) from the rotor into acoustic waves at the 
stator as well as acoustic reflection and transmission properties of both blade rows. In the 
regime of interest for fan design, blade row scattering includes mode and frequency scattering. 
These phenomena have been recognized in the past but never characterized extensively. Hence, 
an important part of the current work has been to explore this scattering behavior and to provide 
a report that can be used to understand the elements of TFaNS and to form a basis for future 
work with more advanced representations of blade row. This report is based strictly on 2D 
analysis and is divided into 2 parts: Part 1 provides theoretical background and a few computed 
examples so that the basic physics can be understood. This part is a reprint of an AIAA paper. 
Part 2 provides an extensive series of scattering curves documenting behavior over a range of 
mode orders and blade passing frequency harmonics. The curves are annotated to aid in physical 
interpretation with reference to theory developed in Part 1. 

The actual code used for scattering in TFaNS has been developed by H. D. Meyer as 
documented in a separate report: “Final Report under Contract NAS3-26618,” December, 1997. 
Meyer treats scattering in a annular duct via a quasi-3D technique: input and output waves are 
3D, blade unsteady loading response is handled by 2D strip analysis, and coupling of the loading 
to duct modes is by the 3D duct Green’s function. The present report and Meyer’s report have 
been coordinated to use the same test cases and a similar notation and presentation format 

The project manager at NASA-Lewis for this work was Dennis Huff. 
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Abstract 

In an earlier paper, a reduced order modeling scheme 
was presented for unsteady acoustic and vortical coupling 
between adjacent blade rows in a turbomachine. The 
essential behavior of the system is governed by modal 
scattering coefficients (i.e. reflection and transmission 
coefficients) of the rotor, stator, inlet, and nozzle, which 
are calculated as if they were connected to non-reflecting 
ducts. The objective of this paper is to identify 
fundamental behavior of these scattering coefficients for 
a better understanding of the role of blade row reflection 
and transmission in noise generation. A 2D flat plate 
unsteady cascade model is used for the analysis with the 
expectation that the general behavior presented herein 
will carry over to models that include more realistic flow 
and geometry. 

It is shown that stators scatter input waves into many 
modes at the same frequency whereas rotors scatter on 
frequency, or harmonic order. Important cases are shown 
where the rotor reflection coefficient is greater than 
unity: a mode at blade passing frequency (BPF) traveling 
from the stator with unit sound power is reflected by the 
rotor with more than unit power at 2xBPF and 3xBPF. 
Analysis is presented to explain this unexpected 
phenomenon. Scattering curves are presented in a format 
chosen for design use and for physical interpretation. To 
aid in interpretation of the curves, formulas are derived 
for special conditions where wavefronts are parallel or 
perpendicular to the rotor or stator chords. 

Introduction 

For the purposes of predicting fan noise, we consider the 
engine to be as sketched in Figure 1 where we identify 
various "acoustic elements": inlet, rotor, stator, and 
nozzle. Much of the fan noise is caused by rotor wakes 
interacting with the stator. Useful prediction models for 
rotor/stator interaction noise have been developed in 
earlier models based on classical analysis methods: 
references 1 and 2, which were early precursors of the 
current model, treated rotor periodic wakes, stator 
loading response, and coupling to duct modes under the 
assumptions of uniform axial flow and no reflections or 
transmission losses from the other acoustic elements. 

^Technical Fellow-Acoustics, Associate Fellow -AIAA 
Copyright © 1997 by D. B. Hanson, published by AIAA with 
permission 
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Figure 1. "Acoustic elements" for a coupled fan noise 
model. 

Later (reference 3), coupling to Eversman’s inlet and 
nozzle radiation codes 4,5 permitted 
computation of harmonic directivity patterns, still 
ignoring reflections back into the source region. 

More recently, emphasis has been on including effects 
of unsteady coupling between the various acoustic 
elements according to the scheme presented in reference 
6. Here, each of the acoustic elements is treated in 
isolation to compute its input/output relationship in terms 
of modal reflection and transmission coefficients. Then 
the elements are coupled at the interface planes by 
solving a linear system that matches the input to each 
element with the output of its neighbors. This coupling 
system is described in detail in the references and is 
summarized in the background section below. Overall 
evaluation of the fan noise prediction system is described 
in a companion paper by Topol (reference 7). The 
purpose of the present paper is to gain a better 
understanding of just one aspect of the coupling physics: 
behavior of the reflection and transmission coefficients 
(or, more generally, scattering coefficients) for the rotor 
and stator blade rows. We go beyond the type of 
information presented in other related studies and 
include scattering into modes and frequencies that are 
different from those of the input waves 

Blade row scattering behavior is computed with flat 
plate cascade theory (an adaptation of the theory of S. N. 
Smith, ref. 8) and steady loading (with the attendant 
turning of the flow) is modeled via actuator disks coupled 
to the blade rows. This simplistic modeling provides a 
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basic understanding of the physics and kinematics of 
scattering and has been used to interpret the behavior of 
the fan noise prediction system. It is hoped that this 
work will also provide a model for future studies with 
more sophisticated computational procedures, i.e. CFD 
methods. This paper can show only a few examples of 
scattering but a supplemental NASA report 9 includes a 
large catalog of scattering curves. Furthermore, a 
separate report by H. D. Meyer 10 includes a quasi-3D 
model with scattering on radial modes. 

The background section below reviews the coupling 
scheme and notation in enough detail that scattering 
coefficients can be defined clearly. That is followed by a 
section that explains the data presentation format and 
reasoning behind it A review of mode kinematics is 
given and then formulas are derived for interpreting the 
scattering curves in terms of wave angles. The M results" 
section presents representative scattering curves with 
discussion and interpretation. Finally, there is an 
analysis of sound power issues, including an explanation 
of how a reflected acoustic power from a rotor can exceed 
the input power. 

Background on Blade Row Coupling and Definition of 
Scattering Coefficients 

This section reviews the coupling scheme used in the 
Theoretical Fan Noise Prediction System (TFaNS) in 
enough detail to define notation used for the scattering 
coefficients that are the subject of this report. Readers 
interested in a complete description of the coupling 
scheme should consult reference 6 or 7. In the coupling 
method, perturbations at the various interface planes in 
the fan are expanded in eigenmode series. The complex 
amplitude coefficients of the series are to be found in the 
coupling solution. Scattering coefficients are defined as 
ratios of these modal pressure amplitude coefficients. 
The scattering curves presented in this paper are in terms 
of power ratios rather than pressure ratios. 

In Figure 1 there are 4 acoustic elements: inlet, rotor, 
stator, and nozzle. These are joined at interface planes 
indexed by r~ 1, 2, and 3. At each of these planes the 
total perturbation in 2D uniform flow can be written as 
the sum of harmonic pressure waves and vortical waves. 
These are indexed, following Smith’s notation, by 7= 1 
(for upstream going pressure waves), 7=2 (for 
downstream going pressure waves), and 7=3 (for 
vortical waves). 

For pressure waves (7= 1,2): 

Pt(x>*>0 = P o £ £a?OU 

n=-oo ir=-«o 

(i) 

Vortical waves (T= 3) axe defined in terms of their axial 
velocity component (rather than directly in terms of the 
vortitity): 


ufOr,#./).* £ ±AKn,k)e* ar « a ' k * x -*fr*'^ 

27 = — OO .£=“<>0 

( 2 ) 

although only the pressure waves are treated in this 
paper. In the above equations, x, cp, and t represent the 
narrow annular duct coordinates and time, n is the BPF 
harmonic index, and k counts the circumferential 
modes. p 0 and a Q are reference pressure and speed of 
sound that can be defined by the user but are usually 
taken to be their sea level, standard day values, a, with 
its indices, represents axial wavenumber and is discussed 

in later sections, xf is the origin for each wave, taken 

at the interface plane where the wave exits an acoustic 
element. (This convention keeps the exponential 
behavior of cut off waves from yielding very large 
scattering coefficients.) Finally, 

IP*,* = + (3) 

is called the kinematic phase since it describes the mode 
spin characteristics. 

m- nB-kV (4) 

is circumferential mode order per the familiar 

Tyler/Sofrin spinning mode theory. 

The significance of equations 1 and 2 is that they 

define the complex amplitude coefficients Af(n, k ) . 

The ordered array of these coefficients is called the state 
vector and its knowledge specifies the perturbations at 
the interface planes. Scattering coefficients are defined 
according to the following notation. 

Af(n',k') <r- Sf.f r (n\ k’\ n, k)Af(n, k) (5) 

This is read: “a wave of type T at interface plane r 
with harmonic index n and circumferential order m = 
nB-kV scatters into a wave of type T at interfa ce 
plane i* with harmonic index n' and circumferential 
mode order m' = n'B-kfV" . The left arrow is used rather 
than an = sign because each input wave scatters into 
more than one output wave. The system equation 
incorporating these scattering coefficients is given in 

figure 2. Bf(n,k)is the source vector, which is 
considered known or prescribed and has the same form 
as the state vector Af(u, k) . The source vector could 

represent wakes from the rotor or it could represent stator 
output waves in an uncoupled environment, i.e. noise as 
calculated for an isolated stator with wake excitation. 
System equation A=SxA + Bin figure 2 is interpreted 
to mean “the state vector A is the sum of the scattered 
waves and the source wave B” 
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System Equation: 
A - S xA +B 


where 

A = State Vector 
B = Source Vector 
S = Coupling Matrix 

Solution: 

A =( 1 -S)" 1 xB 


Figure 2. Generalized coupling scheme for fan noise from ref. 6. 


Certain scattering rules apply. Waves either reflect 
back to the input plane or transmit to adjacent planes, i.e. 

r or r'- r± I. Each wave type T scatters into all 
wave types T'. Scattering on the njc indices is as 
follows: 


inlet: 

n y k 

<- 

n,k 

nozzle: 

n,k 

<r- 

n,k 

rotor 

if ,k 


n^k 

stator: 

n,k ' 

<- 

n,k 


Inlets and nozzles do not scatter on these indices because 
of their assumed axisymmetry. Rotors scatter on 
harmonic order but not on the k index. Stators scatter 
on k but not on frequency. In 3D, there is additional 
scattering or radial mode order. These rules are 
discussed in more detail in reference 6. 


As mentioned in the introduction, scattering 
coefficients are computed by use of Smith’s 2D unsteady 
cascade theory adapted for use with actuator disks that 
address turning of the mean flow. In this application, an 
actuator disk is a plane where the mean flow properties 
jump to represent the gap-averaged change across either 
a rotor or stator. For example, the stator model is built 
up conceptually by starting with the disk separating the 
flow in region 2 (the swirl region) from the flow in 
region 3 (the exit region). The flow on each side of the 
disk is uniform with properties representing the gap- 
averaged flow upstream and downsteam of the stator. 
Mean flow conditions are considered known and are 
prescribed from a separate aerodynamic analysis. The 
acoustic properties of the actuator disk are found by 
requiring that mass and momentum be conserved across 
the disk to first order in the perturbations. (More detail 
is given in ref. 9.) Conservation of mass and 2 
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GtncrqUd wove amplitude 
Incident wove amplitude 


components of momentum across the disk for any wave 
input requires scattering into 3 output wave types (the 
same as described above). Scattering coefficients are 
computed according to equation 5, with scattering only 
on wave type {T into T). Turning of the mean flow is 
associated with the mean loading on the stator. Unsteady 
loading is handled with Smith’s unsteady cascade theory. 
The scattering properties of the cascade are determined, 
again according to equation 3, and the cascade and 
actuator disk are combined into a single acoustic element 
representing the stator by matching the unsteady output 
of each sub-element to the input of the other. 

Format for Scattering Curves and Sample Results 

This section explains the format to be used for reflection 
and transmission (scattering) coefficient curves and 
explains why it is different from formats used in other 
studies. It also presents a sample curve to illustrate the 
major features of interest such as mode scattering, cutoff 
points, and peaks and valleys associated with alignment 
of the waves with the cascade chordlines. Formulas to 
help identify these special features are derived later. 

Smith's Reflection/Transmission Curves - Figure 3 
shows reflection and transmission curves from Smith’s 
report that have often been reproduced by other 
investigators and referenced in reports. The coefficients 
represent ratios of output wave pressure amplitude to 
input wave amplitude for constant geometry, Mach 
number, and reduced frequency while interblade phase 
angle is varied. This is a reasonable format considering 
that the fundamental variables of Smith’s analysis are 


Gap/chord ratio: G = g/c 

Stagger angle: Q 

Mach number: M - W/a 

Reduced frequency: a> = codW 

Interblade phase angle: a 



Figure 3. Smith's reflection and transmission curves. 
G = 1,9 = 60°, A/ = 0.5,(fl =tt /2 


In generation of the curves, the first 4 parameters were 
held constant and the last was varied to vary incident 
wave angle. The frequency chosen by Smith was such 
that only one mode propagates. Hence, the mode 
scattering phenomena highlighted in the present paper 
does not occur, although Smith’s theory includes the 
effect Note from the curves that, in general, the sum of 
the reflection and transmission coefficients is not = 1. 

Rotor/Stator Interaction Noise Analyses - In typical 
rotor/stator interaction analyses, geometry and frequency 
are inter-related in such a way that constant reduced 
frequency is not appropriate for noise studies. Consider 
the geometry and flow sketched in Figure 4 applied to an 
annulus of radius R Often the solidities (or gap/chord 
ratios) of the 2 blade rows are fixed by aerodynamic 
requirements (at least within narrow limits). Also, the 
flow angles and Mach numbers are set by pressure ratio 
and tip speed constraints. In a typical acoustic design 
study, the stator would be varied by changing vane 
number V at constant gap/chord ratio. Since reduced 
frequency depends on chord, it will vary if stator 
gap/chord ratio is held constant. Hence, the strategy for 
the present study was to hold flow angles and Mach 
numbers constant (per the Mach triangle in Figure 4) and 
to vary vane/blade ratio at constant rotor gap/chord G. 
Then reduced frequency can be computed from the other 
parameters using oo = nBCl : 

- one nBClc nJXIR/a c „ nB Mr A, 

<o = — = = = 2n J — /G 

W Ma MR V M ' 

( 7 ) 

Interblade phase angle for the cascade theory is o = - 



Figure 4. Mach numbers and flow angles for rotor and 
stator. Blades & vanes are aligned with flow. 
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2 nnB/V. The role of M T =QR/a is to set frequency. 

The rotor experiences radian frequency co = kVQ . 
In this case the reduced frequency becomes 

coc r kVnc r kVQR/ac r . kV M r f 
r W' M'a M' R B M’ t r 

( 8 ) 

In this case the interblade phase angle is InkV/B , the 
sign difference being caused difference in the sign of 
relative motion of the adjacent blade row. From the 
discussion above, it can be deduced that scattering 
behavior depends only on vane/blade ratio, not on V and 
B individually. 

Figure 5 shows representative curves in the standard 
form of this study. The basic B - V interaction mode 
(i.e. m=nB-kV with n=k=l) was input to the stator from 



Figure 5. Sample reflection curves for stator with 
input of fundmental B-V mode (n,k=l,l). Scattering 
into many modes (n;k f =l;l,2,3.„). 

the rotor side (region 2) and traveled downstream to 
reflect at the stator into many waves (of which 3 are 
shown). The curves give the ratio of power in the output 
waves to power in the input wave with 0 dB representing 
unit reflection ratio. We point out a few features of 
Figure 5 before providing the detailed means to analyze 
scattering curves in the following section. First, note that 
the input wave cuts off at V/B = 2.1 and the k’= 2 and 3 
output waves cut off at V/B = 1.05 and 0.70, 
respectively. Also, note the dip in all of the reflection 
curves around V/B = 0.75. This is the point where the 
input wave is oriented to pass through the stator without 
scattering (wavefronts are normal to the chord lines). 
Finally, the dip at V/B = 1.7 is where the reflected wave 
is oriented for poor coupling to the unsteady loading 
dipoles (again, wavefronts normal to chords). 

Wavenumbers. Wave Angles, and Cutoff 

In this section we develop relationships needed to 
interpret scattering curves in terms of vane/blade ratio 
and the relevant Mach numbers. For example, by 


comparing angles of wavefronts with stagger angles, 
conditions of high transmission or reflection can be 
predicted. Special formulas are derived for the “Venetian 
blind” condition where waves pass through a cascade 
unaffected, the “modal” condition where waves are 
weakly reflected because the loading dipoles are parallel 
to the wave fronts, and the “broadside” condition where 
wavefronts are parallel to the airfoil chords. 

Circumferential Wavenumber and Mode Spin 
Characteristics - In Smith’s original paper, waves are 
specified via complex coefficients multiplying 
exponentials of the general form 

e i(ax+0y+6>t) ^ 

where a and p are axial and tangential wavenumbers 
and co is the radian frequency. For tuibomachinery 
application, the frequency in the stator reference frame is 

co - nBCl (10) 

where Q is the shaft rotational frequency, B is the 
number of rotor blades, and n is harmonic order. Also, 
we wrap the tangential coordinate onto a narrow annulus 
of radius R so that y = R<p and fly = flR <p . For the 
required periodicity in cp , fiR must be an integer. We 

choose 

fiR - -m (11) 

so that the exponential above becomes 

e i{ax-imj>+nBm) ^ 2 ) 

The circumferential mode order m is known from 
standard fan acoustic theory to be 

m=nB-kV (13) 

where V is the vane count and k is an integer index. 
The mode spin behavior can be deduced from the cp, t 
portion of the exponential. If we track the zero phase 
point on a wave, the equation of motion (in the x =0 
plane) is 

m<p Q = nBflt (14) 

from which the mode angular speed = d<p 0 f dt 

is 

^W=— (15) 

m 

Thus, the sign convention chosen above is such that the 
mode spins in the same direction as the rotor shaft for 
positive m (and positive harmonic order n ). Because 
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m can have either sign, modes can be co-rotating (with 
the rotor) or counter-rotating. 

Axial Wavenumber and Cutoff - Smith gives the axial 
wavenumber as 

«u = -r-rh® + ±a V(® + V P ) 2 -(* 2 - u 2 )p 2 

a -U L J 

(16) 

For propagating waves, the 1 subscript and the upper 
sign apply to upstream going waves and the 2 subscript 
and the lower sign apply to downstream going waves. 
Here, U and V are axial and tangential mean flow 
velocity (soon to be replaced by aM x and aM y ). If we 
define a non-dimensional a in parallel with the 
circumferential wave number m, then Equation 16 
becomes 

= — ~~2 x 
l ~ M X 

M x (nBM r - mM y ) ± ^(nBM r - mM y f - (1 - M y )m 2 

(17) 

where M r = QJR/a is the rotor tip Mach number and 
© R/a = nBMr is the non-dimensional frequency 
commonly used in duct acoustics. 

When the square root in Equation 17 is real, the 
wavenumber is real and the mode propagates 
unattenuated in the x direction; when the square root is 
imaginary, the mode decays exponentially. The dividing 
point, where the radicand is zero, is the cutoff condition: 

(oBM Tcutof f - mM y ) 2 =(1- (18) 


The 2 diagonal lines are the cutoff boundaries. The bias 
to the left, associated with a positive M y (swirl in the 
direction of rotor rotation), shows for a given lobe count 
|/n| that backward spinning waves (m <0) cut on at lower 
frequencies than forward spinning waves. Stated 
differently, when the mode spins against the swirl, it cuts 
on at lower frequency. 

Equation 19 can also be written to express the range 
of V/B for cut on modes: 



which is sketched in Figure 7 in a form to aid in 
interpreting the abcissas of the scattering curves. 



Figure 7. Cut on range in terms of V/B ratio and tip 
Mach number. Same abscissa used in scattering 
plots 


The solution for the cutoff Mach number is 


hBMj cutoff — 




m\ +M y 


x m 


which is sketched in Figure 6. 


Wavefront Slopes - In Equation 1 a wavefront is a line 
of constant phase at any fixed time. For example, our + 
P y - 0. Slopes of these lines in the x,y plane are 
(19) = -a //? = 5 x2 f m * From the above expressions, 

this is 


nBM r 




where s m is the sign of m and the upper/lower signs 
apply to up/downstream going waves. 


Figure 6. Mode order - frequency domain showing 
more modes cut on for spin against the swirl. 
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transmitted 



Figure 8. “Venetian blind” condition for stator. Incident 
wavefronts with slope -Mx/My X chord 



Mx 


Figure 10. “Broadside" condition for stator. Incident 
wavefronts // chord. Wavefront slopes = My/Mx. 



wave 

Figure 9. “Modal’ condition for stator: reflected 
wavefront X chord. Wavefront slope = - Mx/My . 


Special Angles and Vane/Blade Ratios 

We expect scattering curves to exhibit certain behavior 
when wavefronts are parallel or perpendicular to the 
airfoils. The most obvious example is when waves 
perpendicular to the vanes impinge on the stator: we 
expect the wave to pass through the blade row without 
scattering as in Figure 8. We call this the “Venetian 
Blind” condition and the figure shows the required slope 
to be -Mx/My. A related, but less obvious, situation 
shown in Figure 9 is called the “modal” condition. Here 
the reflected waves are normal to the chords. With this 
alignment, the unsteady loading dipoles cannot couple to 
the reflected wave and reflection is weak. Again, the 
slope is -Mj/My. 

In the “broadside” condition, waves are parallel to the 
chord lines as in Fig 10. For broadside incident waves, 
we expect high reflection (at least at high frequency). 

The analogous conditions for the rotor are obvious 
extensions of the above. For example, the Venetian blind 
condition is shown in Figure 11. Here, from the Mach 
triangle, the blade chord slopes are -(MT-M y )/M x . 



Figure 11. Venetian blind" condition for rotor. Incident 
wavefronts X chord. Slope = Mx/(Mj~My). 

Now, we set about finding formulas for vane/blade 
ratios that lead to wavefronts either parallel or 
perpendicular to the chordlines. This is done by solving 
Equation 21 for kV/nB in terms of the wavefront slope 
and inserting slopes from the sketches in Figures 8-10. 


kV 

nB 


- = 1 — 


M t 
M y - M X S ± 




( 22 ) 


The sign ±Sm in Equation 21 has been lost because of a 
squaring operation and a new ± has appeared from 
application of the quadratic formula. Hence, we must re- 
evaluate the sign choice based on physical arguments. 

For waves perpendicular to the stator, we showed 
above in Figures 8 and 9 that S= -Mj/My Substitution 
into Equation 22 leads to 
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kV = l _Mr± (23) 

nB M± 1 

To deal with the sign choice, consider the dowstream- 
going wave in Figure 8. The sketch shows that 
downstream-going waves are co-rotating, so that m > 0. 
Since m=nB-kV 9 positive m requires kV < nB so that 
kV/nB < 1. Tins is given by the positive sign in 
equation 23 since the relative Mach number M < 1. 
Thus, the upper/lower sign in Equation 23 applies to the 
down/upsteam going wave. 

For waves parallel to the stator, Figure 10 showed the 
required slope to be M/M x . In this case, substitution 
into Equation 22 gives 
kV 

— = l±M r cosd (24) 

nB 

Here, by similar arguments, the upper/lower sign applies 
to the down/upstream-going wave. 

For the rotor, a more appropriate form of Equation 22 
can be derived: 


nB 


= 1 - 


M t 


kV (M T -My)+M x S±\l+S* 


(25) 


Substitution of slopes from Figure 11 gives the desired 
formulas. All of these results are summarized in Tables 
1 and 2. 


Finally, the cascade response theory can exhibit 
special behavior when the normal distance between 2 
vanes or blades is equal to a half the free space sound 
wavelength. By noting that the normal spacing of the 
stator vanes is gapx cos 0, the frequency is nBQ/2n , 
and the gap is 2nR/V , this condition is easily found and 
is shown in the Tables 1 and 2. 


Table 1. Special Conditions for the Stator* 


Waves Parallel 
to Vanes 

kV 

— = l+A/ r cos# 
nB T 

Waves 

Perpendicular 
to Vanes 

kV _ + M t sin# 
nB~ 1 + M 

y 2 Wavelength 
Normal Vane Gap 

V 

— = 2M r cos 0 
nB 


Table 2. Special Conditions for the Rotor* 


Waves Parallel 
to Blades 

nB 

— = \+M T cos6 r 

Waves 

Perpendicular 
to Blades 

nB ._M r smd r 
kV 1 + M’ 

Vi Wavelength 
Normal Blade Gap 

g 

= 2M t cos#. 

kV T r 

*Upper/lower signs for up/downstream going waves 
Note: by convention, rotor angles are negative 


From these tables, we can find the vane/blade ratios at 
which waves will impinge broadside to the blades (waves 
parallel to blades), will slip through a blade row without 


scattering (waves perpendicular to blades), or will have 
low amplitude when scattered (again, waves 
perpendicular to blades). 


Scattering Results 

The present 2-dimensional study was coordinated with 2 
other studies performed under the LET 4 contract The 
first of these was a parallel study of scattering behavior 
using newly developed quasi-3D capability (ref. 10). That 
study, in turn, used geometry and operating conditions 
chosen to complement evaluation 7 of the TFaNS noise 
prediction system. In the latter case, most of the 
evaluation involved comparison with NASA Lewis’ test 
data from the Pratt & Whitney 22 inch ADP model 
known as “Fan 1”. The quasi-3D study used that fan 
geometry and 2 of the operating conditions, identified as 
the “mid-speed” and “high-speed” cases. The 2D study 
was based on the tip station of the same fan for the same 
2 conditions. These are reported in detail in ref. 9; 
selected results from the high speed case are shown later 
in this paper. 

Fan 1 has 18 blades and 45 vanes. However, in this 
study, vane/blade ratio was varied over the range of 0.5 
to 2.5 while holding the gap/chord ratios constant. The 
geometry and flow were as follows. 


Gap/chord rotor/stator: 
Stagger rotor/stator: 
Pressure ratio: 

Inlet axial Mach #: 


0.962/0.917 
51.3 7 30.5° 
1.278 
0.510 


Rotor rotational Mach #: 0.772 


Given the inlet axial Mach number, the rotor rotational 
Mach number, and the pressure ratio, the remaining 
Mach numbers and flow angles were computed on the 
assumption of zero losses and axial exit flow with results 
shown in Table 3. 


Table 3. 87 5C 

RPM Condition 

Mach 

Numbers 

Upstream 
of Rotor 

Between 

Rows 

Dnstream 
of Stator 

Axial, M x 

0.510 

0.404 

0.388 

Swirl, M r 

0.000 

0.238 

0.000 

Tip, Mr 

0.772 

0.742 

0.733 


Reflection and transmission curves for the stator and 
rotor are shown in Figures 12 and 13. Ordinates are in 
dB expressing the ratio of output power to power in the 
input wave. Thus, 0 dB represents unit reflection or 
transmission. Since emphasis is on rotor/stator coupling, 
input waves chosen for each blade row come from the 
direction of the other blade row. In all cases the input is 
the fundamental interaction mode at blade pacing 
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TRANSMISSION COEFFICIENT - dB REFLECTION COEFFICIENT- dB 


frequency n, k = 1,1, which is to say the m -2? - V 
mode. 

In Figure 12 the stator reflection curve is shown at 
the top and the transmission curve at the bottom. In 
accordance with the scattering rules in equation 6, the 
1,1 input wave scatters into 1,1; 1,2; 1,3, etc. 
Application of Equation 20 to the operating conditions 
on page 8 shows the input wave in region 2 to be cut off 
at V/B = 2.10, which provides the upper limit on the 
reflection curves. In region 3 (downstream of the stator) 
cutoff is at V/B =1.8 as can be seen in the transmission 
curve. The k* = 2 and 3 modes cut off at l A and 1/3 of 
these values, as required by equation 20. At cutoff of the 
1,1 mode, the stator reflects all of the input energy ( 0 
dB ). Application of formulas in Table 1 on page 8 help 
identify prominent features of the curves. For example, 
downstream-going waves are perpendicular to the stator 
chords at V/B = 0.74 according to the middle formula in 
Table 1. This is the “Venetian blind” condition: the top 
curves shows all reflected waves are low in amplitude 
and the bottom curves show transmission of the input 
wave only into itself (no scattering into the other k 9 s). 
The same formula in Table 1 shows the upstream-going 


wave to be perpendicular to the chord at V/B - 1.71. 
This is the “modal” condition indicated by the major dip 
in the reflection curve: the reflected wave is oriented 
normal to the stator loading dipoles and thus cannot 
couple to it. Finally, the bottom formula in Table 1 
shows that a V 2 wavelength fits between 2 vanes at V/B = 
1.28. The dip at V/B = 1.2 in the transmission appears 
to be related to this but the condition is not precise. 

Figure 13 presents the rotor scattering curves in the 
same format. Now the n, k = 1,1 wave at BPF scatters 
into BPF, 2xBPF, and 3xBPF, all at k = 1. From 
formulas in Table 2, a minor dip can be identified at 
V/B = 0.9 where a half wave fits across the blade gap. 
The major dip at V/B - 1.54 is the “modal” condition: 
the fundamental wave input from the stator is reflected 
with wavefronts normal to the rotor chords. More 
interesting, however, is the feet that, with input energy at 
BPF only, the reflected energy is dominated by upper 
harmonics over most of the V/B range. This is 
discussed further in the next section. Finally, the rotor 
transmission curves shown at the bottom of Figure 13 
indicate substantial transmission loss over most of the 
V/B range. However, most of the energy that does get 


STATOR SCATTERING 
INPUT WAVE: *ml, k-1 
OUTPUT WAVES: n«l, l k VI, 2, 3 



ROTOR SCATTERING 
INPUT WAVE: ttml, kml 
OUTPUT WAVES: * 'mJJJ; k ml 





Figure 12. Stator scattering curves 


Figure 13. Rotor scattering curves 
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through the rotor is scattered up from BPF. Modal points 
for the forward scattered waves at 2xBPF and 3xBPF 
are called out at V/B = 0.76 and 1.14. 


S-dLI. M 2 (29) 

a k 


Sound Power Considerations 

In Figure 13 portions of the rotor reflection curves 
pass above the 0 dB level. Since 0 dB represents unit 
power reflection coefficient, this indicates that the rotor 
is reflecting more acoustic power than that in the input 
wave. Values shown at +3 dB indicate that the energy 
reflected at 2xBPF and 3xBPF is twice the input energy 
at BPF. Since there are no simple “conservation of 
energy” principles that apply to moving surfaces, we 
should not be suspicious of this result However, the 
result was unexpected and some analysis is presented in 
this section to explain it 

Power in the acoustic perturbations is evaluated 
according to the energy flux vector given in by 
Goldstein’s Aeroacoustics (ref. 11) 


1 = 


+ u-U 

\Po 


(p 0 u+pU) 


(26) 


This is the instantaneous energy per unit area per unit 
time in the direction given by I. The perturbation 
pressure, density, and velocity are p, p, and u and p G 
and U are the density and velocity of the background 
flow. The energy flux vector is space and time 
dependent. We are interested in its component in the x 
direction. Forming 4 = He and applying p=p/a r 2 lead 
to 


I x =(l+M 2 )pu+-^-p 2 

Pr*r (27) 

+p r a r M x iP' + p r a r M y uv + M x M y pv 


Thus, factor represents the ratio between the modal 

power and modal mean squared pressure. The 
expression for this is derived in the appendix. 

W -MlfnBMrJE 

* ok ~ p T (30) 

( uBM T -mM y ±M x yjE) 2 

where Vi? is the square root that appears in the axial 
wavenumber in equation 17. Since E goes to zero at 
cutoff. Fat does also. In equation 30, the upper/lower 
sign applies to up/downstream going waves. Since we 
dotted I into the unit x vector, we have the necessary 
result that flux on the downstream side is positive and on 
the upstream side is negative (both indicating flux away 
from the source). 

The appendix also derives an expression 
corresponding to equation 30 from the point of view of 
an observer in the rotor coordinate system 

£ F^f (31) 

y r r n k 

where, of course, the pressure coefficients are the 
same in both rotor and stator coordinates. Factor *nk 
is exactly the same as except that nB in the 
numerator of equation 30 is replaced by kV. 

The power ratio plotted in the figures herein is given 
by 10 logio of 


for the power flux in the x direction. Subscripts r in 
equation 27 refer to density and speed of sound in region 
r. When the pressure perturbation given by equation 1, 
and parallel expressions for the associated u and v 
velocity perturbations, are inserted into equation 27, the 
expression can be first averaged over time and then 
averaged over the duct cross sectional area (over q>, in 
this 2D case) to find the average power crossing any x = 
constant plane in the duct. Henceforth, this quantity is 
called “the power” and is in the usual form showing that 
the total power is the sum of the modal powers. 

^ = X F nk\*nkf (28) 

y r T n k 


gg _ F a'k\^D‘k' 

F ak\A akf 


(32) 


which is the power of the scattered wave divided by the 
power of the input wave. An observer in the rotor frame 
would compute a different value 


<H = 


Frik‘\A n 'k’ 

F nk\Aj 


(33) 


The ratio of these ratios is = tikjn k 1 . Since the 
rotor does not scatter on Jc, the ratio is simply 


where y is the ratio of specific heats for air. A r = a/a„ and 
P r = P/Po- Only propagating modes are included in the 
summation. Details of the derivation are given in the 
appendix. A related expression gives the mean squared 
pressure 


5R 


n 


(34) 
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This is to say that the observer in the stator frame will 
compute a reflection (or transmission) coefficient n'/z? 

times what the observer in the rotor frame computes. To 
the rotor observer, the blades are not moving and we 



expect reflection ratios no greater than one. (This has 
been the case in the calculations of this study.) However, 
the change of reference frame results in higher power for 
cases of scattering up in frequency, as indicated by 
equation 34. Note that, in the scattering curves shown 
above, scattering of BPF into BPF did not exceed unit 
power ratio. 

The physical interpretation of the above discussion is 
that the energy in the returning waves is increased by 
motion of the reflecting surfaces. 

Concluding Remarks 

The simple cascade model applied in this paper shows 
how acoustic reflection and transmission of blade rows 
involves scattering into modes and frequencies different 
from those of the input waves. This complex behavior 
can be interpreted with the aid of formulas derived herein 
that help identify peaks and valleys in the scattering 
curves based on orientation of wavefronts relative to the 
scattering blade chords. In particular, the vane/blade 
ratio where wavefronts are perpendicular to the chords 
provides conditions where waves can pass through rotors 
or stators without scattering and also conditions where 
incident waves are not reflected because of dipole 
alignment. It was shown that rotors tend to scatter up in 
BPF order and that reflected energy can exceed incident 
energy by a substantial amount. Although the analysis of 
this paper was based on a 2D flat plate cascade model, 
the kinematic and energy relations are more general and 
will cany over to more accurate representations of the 
cascades and flow. 
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Appendix - Sound Power Formulas 

Equation 27 gives the power per unit area in terms of 
the time dependent pressure and its associated velocity 
components. To derive the expression for the time 
average of this averaged over the duct cross section, we 
insert the series for p from equation 1 and the associated 
series for u and v. Then we average over time, 
integrate over the spatial angle qp, and re-arrange the 
expression to arrive at equation 30 of the main text. 

The series for pressure is 

p(x,*,0 = Po Z (35) 

n k 

and the associated equations for axial and transverse 
velocity perturbations can be written in parallel form 


u(x, I A ak U nJ(e l[a '** + e'*+ +aBQt] (36) 

n k 

Y J A nk V nk e i ^ a ^^ aBat > (37) 
n k 


Here, = -( nB - kV) = -m is the negative of the 

spinning mode order. The axial wavenumber can be 
written 



_ _ 1 | 

M x q±yfE] 

(38) 



in which 





q - dBMj + Myfink 

(39) 

and 

E = g 2 -(1- 

-Mb 

(40) 

Smith’s paper (ref. 8) gives the relationships between 
velocity components and pressure in the pressure waves. 

In our notation, these are 




II 

U-* J 

(41) 

and 

II 

K A-nk ) 

(42) 


where we have defined 

A*nk~ klBMj ^^x a nk + ^yPnk (43) 


In order to apply the derivation to either the rotor- or 
stator-fixed frame, we modify equation 27 slightly by 
adding a prime to M y 

( 1 + )pu+ — p 1 + uv + M x M’ y pv 

Pr*r 

(44) 

Thus, we can set M f y ~ M y , for the observer in stator- 
fixed coordinates, or to M y -M T for the observer in 

rotor-fixed coordinates. Note in equation 44 that, for 
zero Mach number, the intensity becomes I x = pu , the 

familiar result from classical acoustics. 

In equation 44 each term involves a product of 
perturbation quantities, e.g. pu. Since the series 
representations of the perturbations are pure real (by 
virtue of using double sided series in n), we can just as 
well write this as p*u, which is more convenient for the 
ensuing manipulations. When the expressions above for 
p, u, and v are substituted into equation 44, we find 
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Po a o fjr 

r 2 ^ 


S X X X 

i? n' k k* 


xe J l( a nk~ a ii f *') x HfiBk-Po'Jc')# + ( n ~ n, )8& t ] 


(46) 


where y is the ratio of specific heats for air and arose 
from p a al =y p 0 . Also, 


Fno-kk' = (1 + A/j 




+ Af x + A/J g - tt * g ' n '*' 
l ^nk^n'k‘ J 


+M\ 


, ( *nk£ 


n'k' 


V ^■nk^-n'k' J 


+ M x M' y 


-Pn'k' 

*n'k' 


(47) 


When equation 45 is averaged over time, only terms 
with n = n f survive and the <p exponential becomes 
exp i(k - k')V0 . Averaging over cp leaves only terms 


with k-k 1 . Thus, the exponential in equation 45 
disappears completely and the quadruple sum is reduced 
to 


Y r * 


^nk\ 


* n 


where now 

F ak = (\ + M x 

( 


+M'y 


an kPnk \ MM ,{zhL 

#nk ) '{ink 


~T 

l X nk , 

\ 


(48) 


(49) 


After applying the definitions for the a’s, p’s, and X’s 
and considerable manipulation, this can be written as 


Fnk = 

±(1 +M y f} n j c ) + M' y /3 nk ]JE 

{nBMj +Myf$ D ] c i.M x ‘>J~Ey^ 

(50) 


For application to the observer in the stator reference 
frame, we set M' y = M y , in which case 

- ^l-M 2 x fnBM T yf£ 

(nBMj- + MyPrt ± M x -JE) 2 
which is the result quoted in the text in equation 30. For 
the observer in the rotor frame, we set M' = M y - M T 

and we place a hat on F to distinguish the rotor result 


Fnk = 


+i\-M 2 x ) 2 kVM T jE 
(aBM T + M y p nk ± m x 4e ) 2 


(52) 


The upper/lower signs apply to up/downstream going 
waves. On the upstream side, the negative sign indicates 
power flux in the negative x direction: away from the 


source region. The only difference in these expressions 
is the factor nB for the stator where the rotor has kV. 
Thus, observers in the different reference frames will 
differ in their calculation of sound power for any mode by 
the factor nB/kV. 
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PART 2 

COMPENDIUM OF SCATTERING CURVES 
(with Annotation and Discussion) 


Summary 

This part of the report lists the geometry and flow conditions for the scattering study, 
provides tables of special vane/blade ratios keyed to formulas from Part 1 for physical 
interpretation, and presents scattering curves for a range of mode orders and BPF harmonics at 2 
operating conditions. 

Scattering curves are in the format introduced in Part 1: reflection and transmission 
coefficients as power ratios in dB versus vane/blade ratio. The curves are annotated showing 
vane/blade ratios where various input and output waves are either parallel to or perpendicular to 
the chords of the scattering blade rows. Conditions for channel resonances are also identified. 


Conclusions 

A large number of scattering curves was examined to find if qualitative reflection and 
transmission behavior could be predicted based on intuitive expectations regarding wave 
alignment. For example, it might be expected that waves impinging broadside on vanes would 
reflect or scatter strongly and transmit weakly. This was not found to be the case. Similarly, the 
channel resonance condition (where organ pipe modes are set up between adjacent blades or 
vanes) produced nothing remarkable in the curves. 

Regarding the broadside condition, perhaps the results should have not been a surprise. The 
methodology used here is equivalent to that of S. N. Smith (ref. 8) and reproduces the 
reflection/transmission curves in Figure 3. For this figure, the “Venetian blind condition” is at 
60° where it can be seen that the wave transmits perfectly. The broadside condition is at 30° 
where it can be seen that the transmision loss is about 50% (6 dB); the reflection is very weak - 
less than 10% of the input wave amplitude. Thus, the results presented in this report are not in 
conflict with previous knowledge. 

The condition that does produce predictable behavior is that for wavefronts perpendicular to 
chord lines (i.e., where waves propagate parallel to chords). For input waves, this is called the 
“Venetian blind condition”. Here, incoming waves are not reflected but transmit with unit 
transmission coefficient at the input mode order and with no scattering into other modes or 
frequencies. When scattered waves are perpendicular to chords, we call this the “modal 
condition”. Reflected wave energy goes to zero at the modal point and tends to be low in a wide 
range on either side. The result can be a wide valley in the reflection curve when only one mode 
is cut on. 
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2.1 GEOMETRY AND FLOW FOR 2 CASES USED IN STUDY 


A fan with an upstream rotor and downstream stator is represented by flat plate cascades as in 
figure 14. Gap/chord ratio (inverse solidity) matches the tip station of the 22” ADP Fan 1. Flow 


into the rotor and out of the stator is axial; 
flow between blade rows has a tangential 
component associated with pressure ratio and 
tip speed of the rotor. Mean flow conditions 
are specified by inflow Mach number, tip 
rotational Mach, number, and pressure ratio 
based on static conditions of a standard day. 
The remaining flow conditions in the 3 
regions (upstream of rotor, between rotor and 
stator, and downstream of stator) are 
computed under the idealized assumption of 
2D isentropic (lossless) flow, as explained in 
Appendix D. The jumps in the mean flow 
conditions occur at the rotor leading edge and 
stator trailing edge. Stagger angles of blades 
and vanes are computed under the assumption 
that they “weathervane” into the mean flow. 



Figure 14. Geometry and flow of fan stage 

Scattering curves were generated at the 2 

conditions from NASA’s test of ADP Fan 1 shown in Table 4. 


. Table 4. Geometry and Flow Conditions for Scattering Study 



Mid Speed 


Corrected RPM 

7031 RPM 

8750 RPM 

Tip Mach Number 

0.614 

0.772 

Inlet Axial Mach Number 

0.404 

0.510 

Pressure Ratio 

1.180 

1.278 

Gap/Chord: Rotor/Stator 

0.962/0.917 

0.962/0.917 

Stagger Angle: Rotor/Stator 

48.8°, 29.6° 

51.3°, 30.5° 




Static Pressure / Pstandard: Region 1 , 2, 3 

0.894, 1.055, 1.089 

0.837, 1.099, 1.152 

Total Pressure/Inlet Ptot : Region 1, 2, 3 

1.000, 1.180, 1.180 

1.000, 1.278, 1.278 

Speed of Sound/ asTANDARD : Region 1, 2, 3 

0.984, 1.008, 1.012 

0.975, 1.014, 1.020 




Mach Numbers 



Axial: M x i, Mx 2 , M x3 

0.404,0.350,0.341 

0.510, 0.404, 0.388 

Swirl: M v i, M V 2 , M V 3 

0.000, 0.199, 0.000 

0.000, 0.238, 0.000 

Tip Rot’ : Mn, M*n> M*n 

0.614,0.600,0.597 

0.772, 0.742, 0.738 
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2.2 SCATTERING CURVES FOR HIGH SPEED CASE 


Part 1 of this report presented reflection and transmission curves for the high speed condition 
defined in Section 2.1. Results were given only for input of the fundamental wave: n = 1, k = 1 
with scattering into the permitted combinations of n‘ = 1,2,3 and k’ = 1,2,3. This section 
presents curves for more input waves: the 9 combinations of n = 1,2,3; k = 1,2,3 for the high 
speed case. Results for the mid speed case are in Section 2.3. In all cases the input waves 
originate in the region between rotor and stator. 

Presentation of the scattering curves in Sections 2.2 and 2.3 is in pairs of facing pages for 
the same input mode as follows. 

Left Page 


Rotor Scattering 
n ’=1,2, & 3; k=l,2,or 3 


Reflection 

Curves 


Transmission 

Curves 


Vane/Blade Ratio 


Right Page 


Stator Scattering 
n=l,2,or3; k’-l,2,& 3 

1 

Reflection 

Curves 



1 

Transmission 

Curves 

1 

Vane/Blade Ratio 


Scattering Curves for High Speed Case Table 5 gives figure numbers for the scattering 
curves and identifies input and output modes for each. Figures 15a and 15b are repeats of 
Figures 12 and 13 in Part 1. 


Tab 

e5. Figure Numbers - 

High Speed Case 


Rotor 

Stator 

Input 

Mode 

n;k 

Output 

Mode 

n’;k 

Figure # 

Output 
Mode 
n; k’ 

Figure # 

1;1 

1,2,3; 1 

15a 

1; 1,2.3 

15b 

msm 


16a 1 

■HMM 

16b 

1;3 

mmm 

17a 

1; 1,2,3 

17b 

2;1 

■McBM 

18a 

mssm 

18b 

2;2 

1,2,3; 2 

19a 

2; 1,2,3 

19b 

\mxm 


20a 

mxk+*m 

20b 

3;1 

■IHdIM 

21a 

3; 1,2,3 

21b 

■a 

mm jam 

22a 

3; 1,2,3 

22b 

mm 

WEESM 

23a 


23b 


15 

















































Tables 1 and 2 of Part 1 gave formulas for vane/blade ratios at which the various input and 
output waves would be either parallel to or perpendicular to the chords of rotors or stators. The 
first step in study of the scattering curves is to apply these formulas and annotate the curves at 
these special vane/blade ratios. This will aid in interpreting the various peaks and valleys of the 
plots. Tables 6 and 7 below list the special vane/blade ratios for the high-speed condition. 


Table 6. Vane/Blade Ratios for Special Angles Relative to STATOR 


mode order 

Waves 1 Vane Chords 

Waves //Vane Chords 

n, k 


Downgoing 

Upgoing Wave 

Downgoing 

1,1 

1.71 

■ 

0.36 

1.64 

1,2 

■ 


0.18 

0.82 

1,3 

| 

| 

0.12 

0.55 

2,1 

3.42 

1.49 


3.28 

2,2 

1.71 

0.74 


1.64 

2,3 

1.14 

0.50 

^ Hi 

1.09 

3,1 

5.13 


1.08 

4.92 

3,2 

2.56 


0.54 

2.46 

3,3 

1.71 

0.74 

0.36 

1.64 


Table 7. Vane/Blade Ratios for Special Angles Relative to ROTOR 


mode order 

Waves 1 Vane Chords 

Waves // Blade Chords 

n, k 


Downgoing 


Downgoing 

1,1 



1.87 


1,2 



I"'' 1 1 


1,3 

0.13 

0.51 

| 

0.23 

2,1 



3.73 

1.37 

2,2 



1.87 


2,3 



1.24 


3,1 

1.14 

4.63 

5.59 

2.05 

3,2 

€.57 

2.32 

2.79 

1.03 

3,3 

0.38 

1.54 

1.87 

0.68 


These points are marked on Figures 15-23 using the following code names whose meaning will 
become clear in the following discussion. 

BSI Broadside input wave 

VBI Venetian Blind input wave 

MO Modal - scattered wave index is shown 

BSR Broadside reflected wave - scattered wave index is shown 

BST Broadside transmitted wave - scattered wave index is shown 

Rather than discussing each curve individually, we studied the 36 sets of curves in Figures 15-23 
to rind what generalizations could be made. There seems to be as much to say about what 
doesn’t happen as about what does. 

For example, nothing special seems to occur at spinning mode order m = 0. This is the 
plane wave condition computed for nB = kV, which intuition might tell us to avoid for design of 
quiet fans. This intuition is probably based on low speed fan design where any choice other than 
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B-V leads to BPF being cut off. However, for high speed fans, cut on with m = 0 doesn’t 
seem worse than cut on at other spinning mode orders. 

In some simple acoustic problems it can be shown that the sum of the reflection and 
transmission coefficients is equal to 1. Simply stated: what is not transmitted is reflected. A 
quick look at the figures reveals that this is certainly not the case for blade row scattering. The 
sum of the reflection and transmission coefficients can be greater than or less than unity. 

We might expect characteristic behavior for “Broadside” conditions (where wavefronts are 
parallel to the chord lines). For example, the BSI (Broadside input) condition sketched figure 1- 
10 might be expected to result in strong reflection and weak transmission. No such behavior 
was found is study of the curves. Also, for BSR (reflected wave parallel to chord) and BST 
(transmitted wave parallel to chord) there was no tendency toward peaks in the curves. 

We examined the scattering curves to find if there is a tendency toward increased reflection 
at higher frequencies where ray acoustics might be valid. No such trend was found, probably 
because we are dealing with individual waves, not rays. 

We looked for- identifiable behavior when channels between adjacent vanes or blades would 
resonate likft an organ pipe. Appendix C gives the conditions for standing waves in these 
channels and the computed vane/blade ratios are shown in Table 8. 


Table 8. Vane/Blade Ratios for Channel Resonances 


Q 

ROTOR 

STATOR 

k= 1 

k=2 

k = 3 

n=l 

n= 2 

n-3 

l 

0.38 

0.19 

0.13 

2.07 

4.15 

6.22 

2 

0.76 

0.38 

0.25 

1.04 

2.07 

3.11 

3 

1.13 

0.57 I 

0.38 

0.69 

1.38 

2.07 

4 

1.51 

0.76 

0.50 

0.52 

1.04 

1.56 1 


In the table, q is the number of cycles in the standing wave pattern. To avoid clutter, these 
points are not marked on the figures but the reader can verify that nothing special occurs at 
conditions which could produce standing waves. 

The conditions where the formulas of Part 1 reliably predict characteristic behavior are 
those for waves perpendicular to the chord lines. For the VBI (V enetian Blind input) condition 
sketched in figure 8, waves pass through the blade row without scattering (except for weak 
effecis of the actuator disks). Thus, there is little reflection of any waves (see, for example the 
tops of Figures 15b or 18b); on the transmission side, the input wave passes through unchanged 
and there is no scattering into other waves (see bottoms of same figures). The VBI condition 
produces a prominent valley in the reflection curves, in some cases, showing its effect over the 
entire V/B range. 

The modal conditions for transmission and reflection (here, the waves are also perpendicular 
to the chord lines) are reliable predictors of low wave strength although the effect is less 
dramatic than the VBI condition because they affect only one mode at a time. Figure 15b, 
however, shows a case where the modal behavior dominates the plot because it affects the only 
mode that is cut on. 
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Rotor Scattering - High Speed Case 
Input Mode Order n=1; k=2 Output Modes: n -1,2,3; k=2 
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Stator Scattering - High Speed Case 
Input Mode Order: n=1, k=2 Output Modes: n=1; k -1,2,3 
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Figure 16b 
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Input Mode Order. n = 1 ; k = 3 Output Modes, n — 1,2,3, k—3 
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Rotor Scattering - High Speed Case 
Input Mode Order: n=2; k=1 Output Modes: n -1,2,3; k=1 
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Figure 18a 
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Stator Scattering - High Speed Case 
Input Mode Order: n=2, k=1 Output Modes: n=2; k’=1,2,3 
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Figure 18b 
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Rotor Scattering - High Speed Case 
Input Mode Order: n=2; k=3 Output Modes: n -1,2,3; k=3 



Figure 20a 
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Stator Scattering - High Speed Case 
Input Mode Order: n=2, k=3 Output Modes: n=2; k’=1,2,3 
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Figure 20b 


?Q 



Input Mode 


10 T 


0 + 


Reflection 

Coefficient 

dB 



-30 + 


-40 + 


-50 -I 1 1 1 1 h 

0.5 1 


10 T 


0 + 


Transmission , 10 .. 
Coefficient 
dB 

- 20 -- 


-30 + 


-40 + 


•60 -I 1 1 1 1 H 

0.5 1 


Figure 21a 


Scattering - High Speed Case 
: n=3; k=1 Output Modes: n -1,2,3; k=1 



VftX 

i 



Vane /Blade Ratio 


Q n 







Rotor Scattering - High Speed Case 
Input Mode Order: n = 3 ; k = 2 Output Modes, n — 1,2,3, k—2 
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Stator Scattering - High Speed Case 
Input Mode Order: n=3, k=2 Output Modes: n=3;k’-1,2,3 
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2.3 SCATTERING CURVES FOR MID SPEED CASE 


Scattering Curves for Mid Speed Case Curves for the mid speed case are presented in 
this section without discussion since they are similar to those just discussed the high speed case. 
Figure numbers and mode identification are given in Table 9 below and special vane/blade ratios 
are shown in Tables 10 and 1 1 . 


Table 9. Figure Numbers - Mid Speed Case 


1 

Rotor 

Stator 

mnmM 

Output 

Mode 

Figure # 

Output 

Mode 

Figure # 

BBBps 

n’;k 


n; k’ 


i;i 

1,2,3; 1 

24a 

1; 1,2,3 

24b 

i;2 - 

1,2,3; 2 

25a 

1; 1,2,3 

25 b 

i;3 


26a 

mmm 

26b 

2;1 

1,2,3; 1 

27a 

2; 1,2,3 

27b 

mm 

1,2,3; 2 

28a 

2; 1,2,3 

28b 

BS9 

mmm 

29a 


■■ 

WBm 


30a 


WBM 

mm 

1,2,3; 2 

31a 

3; 1,2,3 

31b 

mm 

1,2,3; 3 

32a 


32b 


Table 10. Vane/Blade Ratios for Special Angles Relative to STATOR 


mode order 

Waves 1 Vane Chords 

Waves // Vane Chords 

n, k 

Upgoing Wave 

Downgoing 

Upgoing Wave 


1, 1 


0.79 

0.48 


1,2 


0.40 

0.24 


1,3 

1 

0.26 

0.16 

0.51 

2, 1 

2.99 


■ 

3.04 

2,2 

1.50 


1 

1.52 

2,3 

1.00 

0.53 

0.32 

1.02 

3,1 

4.49 

2.37 

: : ,K wgm 

4.57 

3,2 

2.24 

1.18 


2.28 

3,3 

1.50 

0.79 

0.48 

1.52 


Table 11. Vane/Blade Ratios for Special Angles Relative to ROTOR 


mode order 

Waves 1 Vane Chords 

Waves // Rotor Chords 

n, k 

. Upgoing Wave 

Downgoing 

Upgoing Wave 

iEEZSSHIsSiHHi 

1, 1 

mKMSM 



0.72 

1,2 




0.36 

1,3 

0.17 

■BK9HH 

0.55 

0.24 

2, 1 

1.02 

2.84 

3.31 

1.43 

2,2 

0.51 

1.42 

1.65 


2,3 

0.34 

0.95 

1.10 


3, 1 

1.53 

4.25 

4.96 


3,2 

0.76 

2.13 

2.84 


3,3 

0.51 

1.42 

1.65 

0.72 


37 


































































Reflection 

Coefficient 

dB 


Transmission 

Coefficient 

dB 


Rotor Scattering - Mid Speed Case 
Input Mode Order: n=1; k=1 Output Modes: n'=1,2,3; k=1 
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Stator Scattering - Mid Speed Case 
Input Mode Order: n=1, k=1 Output Modes: n=1; k-1,2,3 




Figure 24b 
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Rotor Scattering - Mid Speed Case 
Input Mode Order: n=1; k=2 Output Modes: n -1,2,3; k=2 




Figure 25a 
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Stator Scattering - Mid Speed Case 
Input Mode Order: n=1, k=2 Output Modes: n=1; k’=1,2,3 
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Figure 25b 
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Rotor Scattering - Mid Speed Case 
Input Mode Order: n=1; k=3 Output Modes: n’=1,2,3; k=3 




Figure 26a 


2 


Stator Scattering - Mid Speed Case 
Input Mode Order: n=1, k=3 Output Modes: n=1; k’=1,2,3 
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Rotor Scattering - Mid Speed Case 
Input Mode Order: n=2; k=2 Output Modes: n -1,2,3; k=2 
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Stator Scattering - Mid Speed Case 
Input Mode Order: n=2, k=2 Output Modes: n=2; k -1,2,3 
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Figure 28b 
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Rotor Scattering - Mid Speed Case 
Input Mode Order: n=2; k=3 Output Modes: n’=1,2,3; k=3 




Figure 29a 
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Stator Scattering - Mid Speed Case 
Input Mode Order: n=2, k=3 Output Modes: n=2; k’=1,2,3 




Figure 29b 
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Rotor Scattering - Mid Speed Case 
Input Mode Order: n=3; k=3 Output Modes: n -1,2,3; k=3 
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Stator Scattering - Mid Speed Case 
Input Mode Order: n=3, k=3 Output Modes: n=3; k -1,2,3 
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APPENDIX B 

ACTUATOR DISK THEORY 


The role of actuator disks in the scattering theory is to account for turning of the mean flow at the 
rotor and stator through conservation of mass and momentum. In applying the conservation 
equations, the flow quantities are written as sums of steady and perturbation (or unsteady) parts. 
When the steady and unsteady parts are separated, the steady flow is considered prescribed (i.e. 
known from a separate aerodynamic analysis. These steady equations provide coefficients for the 
unsteady equations so that jumps in the mean flow cause jumps in the perturbation flow. The 
jumps lead to reflection and transmission coefficients that are found by inverting a linear system. 

There are 3 wave types: 

T=1 : Upstream-going pressure waves 

T=2: Downstream-going pressure waves 

T=3: Vortical (downstream-going) waves 

and 3 conservation equations: 

Conservation of mass 
Conservation of axial momentum 
Conservation of tangential momentum 

These are satisfied on a mode-by-mode basis using the standard wave set from the main text. 

In the following sections, we present the standard wave sets, apply them to the conservation 
equations to develop the linear system, and show the solution method. 


STANDARD WAVE SET 

Forms of the standard waves are given in terms of the defining components and the associated 
components. Pressure waves are defined in terms of pressure; the velocity perturbations are the 
associated components. Vorticity waves are defined in terms of their axial velocity components; 
the tangential velocity is the associated component. 

The defining form for the pressure waves is 


T — 1,2 Pt{**M = Po £ | \At{n,k)e i{aHn ’ kXx - x t )+ '>'°’ k] (B-l) 

J7=— 00 Jc—~CO 

Vorticity waves are defined in terms of their axial velocity component 

T = 3 = f f (B-2) 

D=- aoJc=-co 
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where 


¥ ok = P nk<t> + nBClt 


and 


p^^-inB-kV) 


(B-3) 

(B-4) 


The associated axial and tangential velocity components are given by 


7 = 1,2 




n — — oo k=- oo 


and 


7 = 1,2,3 


Vf (*,*/) = z (B-6) 

D~~PO Jc ——00 


The superscripts r refer to conditions on either side of the actuator disk: r = a on the upstream 
side of the actuator disk and r = b on the downstream side. In the above equations x denotes 
axial distance normalized by radius R of the narrow annulus so that the a’s are also non- 
dimensionalized by R. By applying the momentum equation, or looking up formulas in Smith’s 
paper, the coefficients relating the defining and associated components are easily shown for the 
pressure waves to be 


7 = 1,2 


and for the vorticity waves 


Uf, 


Tnk 


< 


Po 

Pa a o a a 



Po 

Pa a o a a 


T \ 

~ a Tnk 
^ % Tnk ) 


K ^Tnk ) 


We have defined 


7 = 3 


Vfnk = 


~ a Tnk 

Pak 


*Tnk ~ Tip ,r + ^x,r a Tnk + ^s,rPnk 


(B-7) 


0 - 8 ) 

0-9) 


where Mjj p r = (27? / a r is the tip rotational Mach number based on the local sound speed and 
and are the axial and swirl Mach numbers. 


CONSERVATION OF MASS 

J. H. Horlock’s book Actuator Disc Theory (McGraw-Hill, New York, 1978) gives the 
conservation equations applied here. The unlinearized mass flux is pu . Matching this across the 
actuator disk can be expressed 


<J>u) a = ( pu)b 


O*io) 



We write this in terms of steady and unsteady parts (with the tilde’s): 

( Pa Pa X^a a a ) — (P b Pb X^£ ^ b ) 

and focus on the first order unsteady terms, which become 

^ U a ~ __ U b _ 

Pa u a +—Pa=Pb u b +~Pb 


(B-ll) 


a b 


'y 

after applying p = pa . Now, we define 


C = ^ 

Po 


u a 


PaUa+—Pa 


(B-12) 


(B-13) 


for matching at the actuator disk with a similar form C b for region b. The normalization via a</p 0 
is for convenience. a 0 and p 0 are reference values for the sound speed and ambient pressure 
which will drop out of the formulas later. 


We interpret equation 13 to apply on a mode-by-mode basis and express the perturbation 
quantities as the sum over the 3 wave types 


Po 

Po 
a n 


Pa a o 


Po 

Pa a o a a 


~«1 

J 


M x , a 

+ ~~ Po 


Pa a o 


Po 


P a a o a a 


~«2 

y 


+ M **„ 

+ Po 

a* 


+ ~[Pa a o] A 3 

Po 


(B-14) 


where the first line comes from the upstream-going pressure wave, the second from the 
downstream-going pressure wave, and the third from the vorticity wave. Because we are 
matching on a mode-by-mode basis, we have dropped the n,k subscripts and exp iy/ ^ , which 

are common to all terms. All of the matching is at x =0; reference planes are shifted later via the 
exponentials in equations B-l and B-2. We compress the notation and write equation B-14 as 


C a ~ C \a A \ + C 2a A 2 + C 3a A 3 


where 


(B-l 5) 


c _ a o 
C 1 a ~ 


a a 


s~> — a ° 

c 2 a ~ 


a* 


-of 


~«2 

A 


+ Mi 


+Mi 


(B-l 6) 


c - v Ejl 

C 3 a-Y 


( a. ^ 


Po \ a a) 
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•y 

where y , the ratio of specific heats, entered from p Q = yp 0 a 0 . The same argument leads to the 
conserved quantity on the b side of the actuator disk : 

C b -C\b-A\ + C 2b A 2 + ^ 2b A^ (B-17) 

Now we equate C a and Cb . Furthermore, we identify the 3 waves approaching the actuator disk 
as input waves and place them on the right hand side of the equation and the 3 waves leaving the 
actuator disk as scattered waves and place them on the left side. The input waves are the 
upstream-going pressure wave on the b side and the downstream-going pressure wave and the 
vorticity wave on the a side. The scattered waves are on the upstream-going wave on the a 
side and the downstream-going pressure and vorticity waves on the b side. 


C\ a Af -C 2b A 2 -Cy,A 2 - C\bA\ -C 2a A 2 -C^A* (B-18) 

This is the first of the actuator disk equations and represents conservation of mass for the implicit 
n, k mode. 


CONSERVATION OF AXIAL MOMENTUM 

We treat conservation of axial momentum in similar fashion. The unlinearized conservation 
equation is 


Ls\na m + {p + pu 2 ) a = (p + pu 2 ) b (B-19) 

where the first term is the axial component of loading on the disk. The critical step here is to 
separate the total loading on the blade row into the steady part, which is handled via the actuator 
disk, and the unsteady part, which is handled via unsteady cascade theory. Then, when we 
express the variables in equation B-19 in terms of steady and unsteady parts, the unsteady loading 
(on the actuator disk) is by definition equal to zero and the first order perturbation equation is 

(1 + M 1 xa )p a + 2 p a U a u a = (1 + M 2 b )p b + 2 p b U b u b (B-20) 


We define 


F,=—\(\ + Ml a )p a +2p a U,u, 

Pc 1 

and expand it in the 3 wave types, as before 

F a ~ F\A\ + F 2a A 2 + F 2a A* 


(B-21) 


(B-22) 


where 



(B-23) 


( a\ 


^l a =(l + ^i) + 2M ;ra 
F 2a =(\ + M 2 xa ) + 2M xa 


~<*1 

A y 

-«2 ^ 
V A J 


Fla = 2r 


-Pa a o 
/>o 


M 


xa 


The equations for the b side of the actuator disk are the same with b ' s substituded for the a’s. 
By using the scheme described above for identifying source waves and scattered waves, the 
equation for conservation of axial momentum becomes 

F„A{ - F lb A$ - F Jb Al = Fa. 4,* - F 2a A$ - F„A$ (B-24) 


CONSERVATION OF TANGENTIAL MOMENTUM 

The unlinearized equation for transverse (7 direction) momentum is 

-Lcosan + ( puv) a = ( puv) b (B-25) 

Using the same arguments as in the discussion of axial momentum, the first order conserved 
quantity can be written 


G a ~ \Pa G a^a + PaY^a^a xa^ saPa\ 

Po 


or, in terms of the 3 wave types, as 

G a = G \ a A \ + G 2a A 2 + G 3a A 3 


where 


G i a 
G 2 a 

G 3a 


= M 


xa 


= M 


xa 


K A J 
( ZlnL 

\ A ) 


+ M. 


sa 


+ M 


sa 


f 

~ g l 

\A j 

f 

~ a 2 
kA j 


^ xa^sa 


M xa Mga 


— r Fa_^o_ 
Po a a 



( 3 ^ 

“ 


"«3 

+ M 

K Pnk j 

sa 


(B-26) 


(B-27) 


(B-28) 


and the conservation equation for tangential momentum has the same form as the other 2 
conservation equations 


G \ a A l - G 2b A 2 - G 3b A 3 = G \b A \ - G 2a A 2 ~ G 3a A 3 (B-29) 
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SOLUTION OF LINEAR SYSTEM 


By “solution” here, we mean finding the scattered waves for specified input waves. To this end 
we assemble equations B-18, B-24, and B-29 into matrix form 


C\a 

~ C 2b 

~ C lb 

4 


C\b 

G 2a 

-Cla 

1 

^ O'* 

1 

F\a 

~Flb 

~Flb 

4 

= 

Fib 

-F la 

-F 3a 

4 

4la 

~ G lb 

~ G lb. 

4 


G xb 

~G 2 a 

~ G la. 

4 _ 


Solution to this equation is 


, 

i 


C\a 

-Gib 

-Gib 

-1 

' Git 

G la 

Gla 

, 

1 

4 

= 

Fla 

-F 2b 

-Fib 

X 

Fib 

~Fla 

~Fla 

4 

\ 

1 


Pi a 

-Gib 

-Gib. 


Gib 

-Gia 

-Gla. 

4 


or simply 


'4' 



*12 

*13 _ 

A 6 ' 

4 

= 

*21 

*22 

*23 


i 


*31 

*32 

*33. 

1 

iu 

u> ^ 


(B-31) 


(B-32) 


APPLICATION TO ACOUSTIC ELEMENTS 

The objective of this analysis is to create rotor and stator acoustic elements that can be 
represented in the coupling system of Figure 2 of Part 1 . For example, when the stator cascade 
(with uniform flow corresponding to region 2) is combined with an actuator disk at the trailing 
edge to turn the flow from the 0 direction to axial, the result is the stator acoustic element. This 
is done by using equation B-32, which amounts to 3 equations for the input/output behavior of 
the actuator disk, and a similar set of 3 equations for the input output behavior of the stator 
cascade. Setting the output of the cascade (on the downstream side) to the input of the actuator 
disk (on the upstream side) and the output of the actuator disk to the input of cascade permits us 
to eliminate 3 equations. The remaining 3 equations represent the input/output characteristics of 
the combined stator element. 
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APPENDIX C 

CHANNEL RESONANCE CONDITIONS 


When waves at the same frequency travel upstream and downstream in the channel between 2 adjacent 
vanes, a standing wave pattern forms. If the distance between nodes matches the channel length (the 
chord), a resonance could occur that might have an effect on reflection and transmission curves. In this 
appendix, we derive formulas to predict the vane/blade ratios for this condition for stators and rotors. 

From equation 16 in Part 1, we can deduce that the axial wavenumbers for 1 dimensional waves in a 
flow of Mach number M are given by 


, ( D/a 

“'• 2=± I7m 


(C-l) 


Then the sum of a wave travelling upstream and a wave travelling downstream can be written 

p = cos(a jx + cd/) + cos(a 2 x + cor) (C-2) 

Since we are only finding the standing wave condition and not attempting to be general, we use unit 
amplitude waves. The mean square of the pressure can be obtained by squaring this expression and 
averaging over time with the result 


p 2 =l + cos 


2 — 
S_ 

l-M 2 

J 


(C-3) 


The resonance condition occurs when the phase of the cosine equals 2 ji, 4k, etc. at x = b , i.e. for 

_ a)b 
2 — 

a T = 2n 9 

l-M 2 (C-4) 

where q is an integer. For the stator the radian frequency is d) = nBfl and the chord can be 
expressed as 

1 M 

~G V (C-5) 

where G is the stator gap/chord ratio. Substitution into the Equation C-4 gives the desired formula for 
the stator 


V Mr 

nB (1 


Similar analysis for the rotor leads to 


B 


M T 


kV (1 -M /2 )|G r 


(C-6) 


(C-7) 


where M’ is relative Mach number in the rotor blade passages and G r is rotor gap/chord. These 
formulas are used in Section 2.2 of the text to generate Table 8. 

62 



APPENDIX D 

CALCULATION OF MEAN FLOW PARAMETERS 


Analysis in this report recognizes 3 regions of mean flow: 

Region 1 Upstream of rotor axial 

Region 2 Between rotor and stator swirling 

Region 3 Downstream of stator axial 

Background (or mean) flow values of Mach number, pressure, and speed of sound in each of these 
regions are required for the acoustic and unsteady aerodynamic models. This appendix presents the 
scheme for calculation of the mean flow, including the FORTRAN code listing and a sample case. 
Much of the documentation is included with the listing. 

Flow regions and some of the associated notation is sketched below. 


Rotor Stator 



Region 1 
upstream 


Region 2 Region 3 

mid-stage downstream 


M x i M x2.My2 M x 3 

Pi P 2 P 3 

T l T 2 T 3 


Operation of the system is as follows. Upstream total pressure and temperature are fixed at sea level, 
standard day values. The user specifies inflow axial Mach number, rotor rotational Mach number 
(based on speed of sound from region 1), and pressure ratio. The remaining quantities are then 
computed with ID isentropic (lossless) flow equations and the standard equation for action of the rotor: 

C 

A(tangential velocity) 2 x A(total temperature) 

wheel speed 


where Cp is the specific heat at constant pressure. Steps in the procedure can be followed easily in the 
code listing. 
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Calling Routine 

The code is written as a subroutine. The calling main program immediately below was written for 
testing purposes. It is followed by the background flow subroutine. Sample I/O is shown at the end of 
this appendix. 


PROGRAM CALLMEANFLOW 

IMPLICIT DOUBLE PRECISION (A-H, O-Z) 
DIMENSION A (3) , P(3) 

DOUBLE PRECISION MX (3 ) ,MY (3 ) ,MT (3 ) ,M1 ,MT1 


1 

2 

3 

4 

5 


WRITE (*,*) 'Input Ml, MT1 , 
READ(*,*) Ml , MT1 , PR 
CALL MEANFLOW (Ml , MT1 , PR , 
WRITE ( * , 1 ) MX(1) , MX (2 ) 

MY (1) 

MT(1) 

A(l) , 

P(l), 


WRITE (* , 2) 

WRITE (* , 3 ) 

WRITE (* , 4) 

WRITE ( * , 5 ) 
FORMAT ( IX , 1 MX1, 
FORMAT (IX, 'MYl, 
FORMAT (IX, 'MT1 , 
FORMAT (IX, ' Al, 
FORMAT (IX, ' PI, 


MY (2) , 
MT(2) , 
A (2 ) , 

P (2) , 
MX2, MX3 
MY2 , MY3 
MT2 , MT3 
A2 , A3 
P2 , P3 


PR' 

P, A,MX,MY,MT 
MX (3 ) 

MY (3 ) 

MT (3 ) 

A(3) 

P(3) 

', 3F10.3) 

', 3F10.3) 

', 3F10.3) 
', 3F10.3) 

• , 3F10.3) 


) 


END 


Background Flow Code 


SUBROUTINE MEANFLOW (Ml , MT1 , PR , P, A,MX,MY,MT ) 
c.. This routine computes the mean (or background) flow parameters needed for 
c. acoustic calculations in 3 regions: l-> upstream of rotor, 2-> between 
c.. rotor and stator, 3-> downstream of stator. 

c.. Calculations based on rotor action being given by delta T0=U*Cy2/Cp, 
c. . one dimensional isentropic flow, equation of state (to get speed of sound) , 
c. . continuity, and assumption of zero losses. Coding by D. Hanson, Consulting on 
c.. aero theory by M. Devries 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


Regions : 1 


axial 

flow 

upstream 


rotor 

\ 

\ 2 
\ 


stator 

/ 

/ 

/ 


3 


> 


/ 

/ 

/swirl 
/ region 

>/ 


axial 

flow 

downstream 


SET: T01 , P01 , Ml 

as inflow conditions 

SPECIFY: Pressure ratio (PR) and rotor rotational Mach number (MT1) 
COMPUTE: 

Tl, PI, Al P2, T2 P3 , T3 , 


64 



IMPLICIT DOUBLE PRECISION (A-H, O-Z) 

DIMENSION A (3) , P(3) 

DOUBLE PRECISION MX (3 ) , MY (3 ) ,MT (3 ) ,M1 ,M2 ,M3 ,MT1 ,M20LD,M30LD 


c.. Input based to be based on Sea level standard day conditions 


GAM = 1.4D0 
R = 1716. 59D0 
Cp = 6008. 09D0 


Ratio of specific heats 
Gas Constant 

Specific heat @ constant pressure 


Tstd = 518.67D0 

Pstd = 2116. 22D0 

Astd = SQRT (GAM*R*Tstd) 


! Sea level Standard Day - Degrees Rankine 
! Sea level standard Day - pounds/sq.ft. 

! ft/sec (from equation of state) 


T01 = Tstd ! Inflow based on Sea Level Standard Day 

P01 = Pstd ! Inflow based on Sea Level Standard Day 


c.. Compute parameters for Region 1 (upstream of rotor - axial flow) 

FI = 1 + (GAM-1D0) /2 * Ml**2 ! ID isentropic flow 

T1 = T01 /FI ! ID isentropic flow 

PI = P01 / FI** (GAM/ (GAM- 1D0) ) ! ID isentropic flow 

RH01= P1/(R*T1) ! Gas law 

A1 = SQRT (GAM*R*T1 ) ! equation of state 

U = MT1 * A1 ! rotor tip Mach number based on region 1 

Cxi = Ml * A1 ! axial velocity (needed later for continuity) 


c.. Compute parameters for Region 2 (between rotor and stator - swirl flow) 
P02 = PR * P01 ! Pressure ratio from input 

DELTATO = T01 * ( PR** ( (GAM- 1D0) /GAM) - 1D0) ! assumes 100% efficiency 

T02 = T01 + DELTATO 

Cy2 = (Cp/U) * DELTATO ! Equation for action of rotor 


c. . 


10 


Iterate on Mach number and continuity 
M2 = Ml * SQRT (Cxl**2+Cy2**2) /Cxi ! First guess to start iteration loop 
write (*,*) 'Iteration on Mach number in region 2' 


DO I = 1, 10 

F2 = 1 + (GAM-1D0) /2 * M2**2 
T2 = TO 2 / F2 

P2 = P02 / F2** (GAM/ (GAM-1D0) ) 
RH02=P2/ (R*T2 ) 

A2 = SQRT (GAM*R*T2 ) 

Cx2 = Cxi * RH01/RH02 
M2old = M2 

M2 = SQRT (Cx2**2+Cy2 
write(*,*) i, m2 
IF (ABS (M2-M20LD) 

ENDDO 
WRITE (* 

CONTINUE 


! ID isentropic flow 
! ID isentropic flow 
! ID isentropic flow 
! Gas law 

! equation of state 
! continuity 

! save for convergence check 
! new estimate of Mach number 

10 

CONVERGED IN MEANFLOW ROUTINE' 


**2) / A2 

.LT. . 0001D0) GOTO 
*) 'WARNING, M2 LOOP NOT 


write (*,*) • • ! write swirl angle (from axis of rotation) 

THETA = ATAN ( Cy 2 / Cx2 ) 

write(*,*) 'Flow angle = ', theta*180d0/3 . 14159265D0 
write(*,*) ' ' 

c.. Compute parameters for Region 3 (downstream of stator - axial flow) 

T03 = T02 ! Assume zero losses 

P03 = p02 ! Assume zero losses 


c. . 


Iterate on Mach number and continuity 
M3 = M2 * COS (THETA) ! First guess to start iteration loop 

write (*,*) 'Iteration on Mach number in region 3' 


DO I = 1, 10 

F3 = 1 + (GAM- 1D0) /2 * M3**2 
T3 = T03 / F3 

P3 = P03 / F3** (GAM/ (GAM-1D0) ) 
RH03=P3/ (R*T3) 


! ID isentropic flow 
! ID isentropic flow 
! ID isentropic flow 
! Gas law 
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20 
c. . 


A3 = SQRT (GAM*R*T3 ) 

Cx3 = Cx2 * RH02/RH03 

M3old = M3 

M3 = Cx3 / A3 

write (*,*) i, m3 

IF (ABS (M3-M30LD) .LT. .0001D0) 


! equation of state 
! continuity 

! save for convergence check 
! new estimate of Mach number 

GOTO 20 


ENDDO 

WRITE (*,*) 'WARNING, M3 LOOP NOT CONVERGED IN MEANFLOW ROUTINE' 
CONTINUE 

Prepare values to return to calling routine 
A ( 1 ) = Al/Astd 
A (2) = A2/Astd 
A (3) = A3 /Astd 
P(l) — Pl/Pstd 
P (2 ) = P2/Pstd 
P (3 ) = P3/Pstd 


MX (1) 

= Ml 


MX (2 ) 

= M2 * 

COS (THETA) 

MX (3 ) 

= M3 


MY (1) 

= ODO 


MY (2 ) 

= M2 * 

SIN (THETA) 

MY (3 ) 

= ODO 


MT(1) 

= MT1 



MT (2 ) = MT1*A1/A2 
MT (3 ) = MT1*A1/A3 


RETURN 

END 


Sample Case 


Type a.out to get: 

Input Ml, MT1 , PR 

Enter values for axial Mach number, tip rotational Mach number, and pressure ratio: 

.510 .772 1.278 


Hit “Enter” to do the calculation with the result: 


Iteration on Mach number in region 2 

1 0.49069069566595 

2 0.47318442132157 

3 0.46959417475541 

4 0.46887728952648 

5 0.46873493006882 

6 0.46870669130413 

Flow angle = 30.501298137168 

Iteration on Mach number in region 3 

1 0.39084584049871 

2 0:38850434236162 

3 0.38809172279216 

4 0.38801929795917 


MX1, 

MX2, 

MX3 

- 

0.510 

0.404 

0.388 

MY1, 

MY2 ; 

MY3 

= 

0.000 

0.238 

0.000 

MT1, 

MT2 , 

MT3 

’ = 

0.772 

0.743 

0.738 

Al, 

A2, 

A3 


0.975 

1.014 

1.020 

PI, 

P2 , 

P3 

= 

0.337 

1.099 

1.152 
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